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ABSTRACT 


Ray tracing of the Jovian magnetosphere in the low frequency range (1--40 
MHz) has resulted in a new understanding of the source mechanism for lo depen- 
dent decametric radiation (DAM). Our three-dimensional ray-tracing computer 
code has provided model DAM arcs at 10® intervals of lo longitude source posi- 
tions for the full 360“ of Jovian system III longitude. In addition, particu- 
larly interesting arcs have been singled out for detailed study and modelling. 
The results of this study indicate that DAM arcs may be categorized according 
to curvature — the higher curvature arcs are apparently due to wave stimulation 
at a non-constant wave normal angle, (j** The 4^{f) relationship has a signature 
that may be common to most of the higher curvature arcSo The low-curvature 
arcs, on the other hand, seem to be adequately modelled with a constant wave 
normal angle of close to 90“ . These results imply that for higher-curvature 
arcs observed far from Jupiter (to diminish spacecraft motion effects) the 
electrons providing the gyro-emission may be relativis tically beamed. 
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IN'rRODUCTION 


One of the moat characteristic features of the Jovian high frequency 
radio emissions is the presence of the arc-Xike structures in frequency-time 
spectrograms. These emissions have been observed in the decameter (DAM) fre- 
quency range of 1 MHs to 40 MHz and the arcs occur as single-valued and 
double-valued in frequency and as a series of closely-spaced (nested arcs in 
time (see Figs. 1 and 2). Careful study has disclosed the arc-like structures 
even on ground-based spectrograms (Boischot et al. » 1981). 

The radiation is observed both in the right and left hand optical polari- 
zations as if coming from source regions in the northern and southern heri- 
spheres. Boischot and Lecacheaux (1981) have noted a remarkable repeatability 
of features at certain lo phases. The arcs can be described as being vertex- 
early or vertex-late if the "nose frequency" occurs at earlier or later Jovian 

longitudes respectively. In addition the arcs are labelled as "great arcs" if 
their frequency extent is ~5 MHz to ^20 MHz . Sources are described as being 

either A, B, C, or D if arcs are observed in the following Jovian system III 

longitude ranges: source A or main source (190° < CML < 290®); source B or 

early source (60° < CML < 190°); source C or late source (290° < CML < 30°); 

source D (30° < CML < 60°). 

A study of the Voyager I and II data shows that many of the nested arcs 
are not double-valued arcs (double-valued in frequency), but rather a series 
of small-radius-of-curvature arcs. Total time duration of an individual shal- 
low arc typically ranges from three to thirty minutes. The time interval be- 
tween nested arcs can be quite small, sometimes less than one minute. The 
hi gher-radius-of -curvature arcs are longer in time duration (possibly >2 hrs) 
and sometimes are nested in groups as well. 
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Various source mechunisms have been proposed to explain decametrlc arcs. 
Many Investigators attribute the emissions to precipitating electrons along 
the lo flux tube which then Interact with the Jovian Ionosphere and upper at- 
mospheric plasma via a gyro-emlsslon mechanism (cf Neubauer et al, 1980). 

The resulting radiation Is envisioned as being emitted In a hollow cone 
of half angle, <J*- with respect to the B-fleld. Warwick et al. , (1981) suggest 
the arc structure Is due to anomalous B-fleld Islands at the base of the lo 
flux tube. Bolschot (1981) suggests an Interesting Interference effect with 
the lo torus as the source of the arcs. In this paper we suggest propagation 
effects through the Jovian Ionosphere and magnetosphere as the arc-producing 
mechanism. 

Some preliminary work has been done to describe the growth rate of the 

t 

earth's AKR emissions, analogous emission to the Jovian DAM radiation (cf. 
Omidi and Gurnett, 1982; Wu 1982; Calvert, 1982). These authors' find- 

ings suggest that for conditions as might be expected at the foot of the lo 
flux tube (i.e. <j^/a^<1) emission should occur at large (^60°) angles with 

respect to the magnetic field (large wave normal angles). 

In this paper three-dimensional ray-tracing of the Jovian DAM radiation . 

is performed using a realistic Jovian plasma model (Sentman and Goertz, 1979), 

and the 0-4 magnetic field model (Acuna and Ness, 1979). The purpose of this 

paper is to perform three-dimensional ray-tracing of DAM radiation from source 

positions along an lo flux tube at frequencies close to the right-hand 
extraordinary (RX) mode cutoff frequency, f^x = fg/2 + fg^+ fp^ (fg = gyro- 

frequency; fp = plasma frequency) . Jovian arcs are produced with minimal 

assumptions about the source mechanism. We suggest it is possible to explain 
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the low-curvature arcs as due to the mechanism of multiple reflections of 
Alfven waves as presented by Neubauer (1979) and Goertz (1980), and elaborated 
by Gurnett and Goertz vM981). A constant wave normal angle at each frequency- 
dependent source point sufficiently reproduces the observed results. We ex- 
plain the large sometimes isolated high-curvature arcs as due to source points 
with wave normal angles which vary as a function of frequency in a repeatable 
fashion, and we suggest these latter arcs are a result of more energetic field- 
aligned particles than are the low-curvature arcs. 

B-field and Plasma Models. 

The raypaths calculated by the procedure described above are dependent on 
the magnetospheric model chosen. The description of the Jovian magnetic field 
utilized in these calculations (the 0-4 model) is that published by Acuna and 
Ness (1976); the background plasma, Sentman and Goertz (1978); and the lo 
torus, Warwick, et al (1979). Figure 3 is an abbreviation of this combined 
"Jovian magnetosphere" showing contours of plasma frequency in the 200° -20° 
meridian plane. Suspected source regions in the northern and southern hemi- 
spheres are shaded. Spherical harmonic expansion coefficients are used to cal- 
culate the magnitude of the field Bq, the components Bj-, Bq, and B(() and their 
derivatives. The plasma density along the raypath, and the spatial gradients 
are calculated from empirical fits. 
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Emlsalon ’Mechanism 

In this study the following assumptions have been made concerning the 
Jovian emission mechanism: 

1. Tlie radiation is in the right-hand polarized extraordinary mode. 

2. The source region is at the foot of the lo flux tube and the source 
mechanism is Doppler-shifted gyroemission, 

3. The emission cone is hollow. 

Assumption No, 1 is stongly supported by the data set (cf. Warwick et, 
al.). The strongest emissions are RX, and when left hand polarization is pre- 
sent, the source is most likely the southern hemisphere. For this work the 
emission mechanism assumed is that first presented by Neubauer, (1979); 

Goertz, (1980); Gurnett and Goertz, (1980). The mechanism depicts lo inter- 
acting with the Jovian magnetosphere and generating large amplitude standing 
Alfven waves which propagate northward and southward along the magnetic field 
lines from lo. Each Alfven wa.ve is associated with a pair of oppositely di- 
rected currents, in this model, which flow along the Alfven precipitation path 
between lo and the Jovian ionosphere. Assumption No, 2 is consistent with 
this picture. The degree of Doppler shifting is one of the subjects of this 
investigation. Finally, a hollow emission cone for the DAM is strongly sug- 
gested by the data (Carr and Gulkis, 1969), Analogously to terrestrial auro- 
ral kilometric radiation (AKR) the initial DAM radiation at a given frequency 
is assumed approximately perpendicular to the magnetic field, i.e., the ini- 
tial wave normal angle, (|i, is aproximately 90° (this assumption is later re- 
laxed) . Those rays emanating from a lower Jovian source altitude eire closer 
to the R-X cutoff frequency, and thus are refracted more. While AKF. emissions 
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have filled emiasion cones, Jovian DAM emissions, whose source points are in 
regions of much lower density than AKR emissions, are believed to have hollow 
emission cones. Another important point is that the wave normal angle is a 
function of frequency. In general the higher frequencies with source points 
close to the surface have smaller wave normal angles than do lower frequencies 
with higher altitude source locations, 

RESULTS 

To calcv’late propagation ray paths in the Jovian magnetosphere we have 
used a three-dimensional ray-tracing program. The computer code is based on 
the Stix cold-plasma formulation of the index of refraction and Haselgrove's 
set of first order differential equations amenable to numerical solution. 
Chosen as inputs are: a source position (in radial distance, longitude, and 
latitude) from which waves at a specific frequency, f, are to be launched at a 
specific wave normal angle (WNA), (|j, with respect to the magnetic field. Dif- 
ferent frequencies originating from source points on the same field line or on 
adjacent field lines would generate a set of nested emission cones. Depending 
on the geometry of these cones, imposed by the plasma and by the B-field, the 
spacecraft would detect frequencies v;hich decrease or increase with time as 
these cones sweep over the detector. 

At the source region, the program calculates an index of refraction sur- 
face for R-X mode waves, based on B-field and on background plasma parameters. 
Next, the program takes an incremental step in the direction perpendicular to 
the index of refraction surface, i.e,, in the direction of the group velocity 
or energy flow, and then determines the coordinates of this new point on the 
raypath. Then, another index of refraction surface is calculated through this 
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new point and the steps are repeated. The direction of each ray changes ac- 
cording to Snell's law as it travels through the magnetized plasma. To gene- 
rate an emission cone, thirty-six rays are launched, one every 10® around B at 
the specified wave normal angle and frequency. 

Radiation Sheets of Varying Frequency. 

Figures 4 and 5 present mercator projections of certain radial shells 
centered on Jupiter. The intersections of the calculated, model -dependent 
emission cones with these radial shells are labeled according to frequency. 

In Figure 4, the radial shell is at 6Rj and the lo torus is shown in shaded 
cross-section. Jovian latitude is indicated ',>n the vertical axis and longi- 
tude System III '65, on the horizontal axis. Each point on the figure is the 
intersection of a raypath with the 6Rj shell. In three dimensions, not shown 
here, the rays could be visualized as the bent ribs of an umbrella, while the 
fabric would correspond to the warped emission cone. All of the rays in Fig- 
ure 4 were launched at (J> = 90®. Emission cones of two frequencies, 2 MHz and 
20 MHz, originate from source points in the northern hemisphere labeled "X", 
and emission cones of 2 MHz and 10 MHz radiation originate from source points 
in the southern hemisphere. At 6Rj parts of the emission cones are penetrat- 
ing the model lo plasma torus with no noticeable effect. In fact, the torus 
was "invisible" to all DAM frequencies of this study. From close scrutiny of 
Figure 4, one notices the source points for different frequencies are dis- 
placed in longitude (horizontally), as well as in latitude. Since different 
frequencies are launched near the local gyrofrequency, one expects the source 
points to have different latitudes. The longitudinal separation of the source 
points for two frequencies is an indication of the cizimuthal twist in the 
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field line threading lo. This field lino is shown dashed, as it is projected 
onto the 6Rj shell. Also, notice that the emission cones of differing fre- 
quencies can intersect, and cross over one another. Finally, notice how the 
emission cones directed toward the southern 'fj^misphere are tilted around the 
source point and appear to be rotated with respect to the northern hemisphere 
emission cones. 

The trajectory of the raypaths out to a distance of 100 Rj is illustrated 
in Figure 5, I’he points represent the same rays seen in Figure 4, but now 
traced to lOORj. Comparing Figure 4 with Figure 5, it is easy to see that the 
emission cones have spread out. Near longitude 230®, the radiation from the 
northern hemisphere sources reaches a southern latitude of -60®, while the ra- 
diation emanating from the southern sources reaches a northern latitude of 
45®. The radiation sheets at lOORj intersect the equatorial plane at dif- 
ferent points (more widely separated) compared with the intersections of the 
same sheets with the equatorial plane at 6Rj. Notice that the radiation 
sheets at 100 Rj are smooth and have the same gentle curvature as those at 6 
Rj. Since a portion of the emission cones traversed the lo plasma torus at 
6Rj, one might expect "kinks" to occur in the cones, beyond the torus. No 
evidence of refraction effects, caused by tlie model torus, are seen in the 
lower panel. We conclude that the torus has negligible effect even at 2 MHz, 
the lowest frequency traced, and leads us to surmise that the major propaga- 
tion effects occur near the source regions. 

Closer examination of the projection at lOORj reveals "crossings" of the 
emission cones of different frequencies . This should not be surprising since 
different frequencies are launched from different source points, assuming ra- 
diation occurs near the local electron gyrof requency . The twist in the field 
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lines from one source point to the next produces intersectiows of the emission 
cones of different frequencies. The trajectory of Voyager 1 was slightly 
above the Jovian equator. As the magnetic field rotates, the emission cones 
sweep from left to right over the spacecraft (or, in the Jovian frame, the 
spacecraft flies from right to left, slightly above the equa;i;or) , As the 
spacecraft intercepts these emission cones, the PRA experiment records fre- 
quencies which decrease or increase with time, depending on the nature of 
these intersections, i.e,, the geometry of the emission cones. 

Radiation Sheets of Varying Wave Normal Angle 

The top panel of Figure 6 contains emission cones for 2 MHz and WNA's of 
60», 80«», 100*, and 120°. Rays with ittitial WNA >90° (4<>90°) are initially 
directed toward the planet. This mercator projection was taken at the lOORj 
radial shell. Note that the emission cones "open up" as the WNA increases 
from 60° to 100°, then the radiation emitted at (J» = 120° refracts upwards away 
from the R-X cutoff. The bottom panel illustrates the same effect for emis- 
sion cones of 25 MHz. WNA is one parameter which shifts the intersection of 
the emission cones at the equator (and at the spacecraft) toward earlier or 
later longitudes. 

Radiation Sheets Emanating from Varying Source Longitudes 

Figure 7 shows emission cones of frequency 25 I4Hz, WNA = 90° at the lOORj 
shell. The projected source positions are located at the top of the frame. 

If adjacent field lines are sources of emission at a given WNA, a family of 
"nested emission cones" of that particular frequency is produced. The refrac- 
ting of the emission cones from left to right across the frame may be caused 
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by the azimuthal aoeymmatry i,n the field, which is an Integral part of this 
model. As the Jovian magnetic field rotates, the spacecraft would be illumi- 
nated by this frequency as long’ as these emission cones intersect the space- 
craft position. The spacecraft could "fly out of" this emission when sources 
are absent on adjacent field lines or when the geometry of the emission cones 
produces no Illumination of the spacecraft. 

MODEL ARC RESULTS 

A typical signature of arc quite common in the data is that of the short 
time duration, low-curvature nested arcs which can extend over a frequency 
range of f >30 MHz. (Pig. 2). In order to construct model low curvature arcs 
the dynamics of the problem were considered: Jupiter was rotating; lo was 

revolving; and the spacecraft was moving. To minimize the effects of space- 
craft motion the ray path emission cone projections were examined at large 
radial distance (150 Rj) . During the time of observation of an arc the central 
meridian longitude of Jupiter changes. For each change in lo system III 
longitude of 10®, the spacecraft system III longitude varies by about 13® 

(when the spacecraft is far from Jupiter). In this study we assume that the 
lo flux tube is the radiation source. Initially we consider only the case of 
emission from a series of constant lo longitudes, each successive lo longitude 
10° different from the previous. We also initially assume emission at only 
one wave normal angle, » 90®, Figure 8 is a plot of frequency versus 
spacecraft system III (1965) longitude. On the right in this figure curves 
for three different lo flux tubes are shown: lo longitude = 180°, 190®, and 

200°. Each point on each model curve represents the spacecraft longitude at 
which the chosen frequency intersects the spacecraft latitude at the time of 
observation. Source points for the seven different frequencies chosen are 
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located along each lo flux tube at pooltiono juot above the R-X cutoff fro~ 
quency* For all cases where the wave normal angle, (|/, is approximately 90®, 
which was the case for the low curvature nested arcs modelled, the Doppler 
shift was assumed minimal: f/fg 1.005. This value is just large enough to 
prevent singularities at the source point in the computer code. The series of 
model low curvature nested arcs shown (vertex>-late in this case) qualitatively 
matches a series of nested arcs observed by the FRA experiment on Voyager 2 
when the spacecraft longitude and lo's longitude were dose to those used in 
the model case . One of tlie actual arcs is plotted for comparison in Fig. 8. 
There is a slight discrepancy between the model arcs and the data at the 
higher frequencies. This discrepancy may be due to emission from slightly 
smaller wave normal angles at the higher frequencies. For the case of most low 
curvature nes^ad arcs, the time duration of an arc is so small (several min.) 
that Zv fliiX tube longitude of the source can be assumed constant for the 
duration of the arc. Also shown in Figure 8 are two model vertex-early-nested 
arcs (quite small curvature in this case) and the observed Voyager 2 arc at 
the approximate same spacecraft system III longitude and appropriate lo flux 
tube source longitude. The two model arcs shown bracket the observed arc. 

The model shows good qualitative agreement with the actual data, but the model 
lo longitude and actual lo longitude during the time of the observed arc dif- 
fer by approximately 25® . This discrepancy is probably due to the difference 
between the magnetic field model (0-4) and the actual magnetic field at low 
altitudes near the source point. 

We have modelled arcs for source longitudes at 10® intervals for the full 
360® of System III longitude. The model arcs are all for an assumed Doppler 
shift frequency ratio of f/fg = 1.005 and (|> » 90®. 


To construct the arcs it 
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was assumed that the spacecraft was at 150 Rj and at a latitude of 3»2®« 

Again, each point; on each model curve represents the spacecraft longitudo at 
which the chosen frequency intersects the approximate spacecraft latitude at 
the time of observation. Hence each point was determined by intersecting the 
approximate spacecraft latitude with the model-generated emission cone surface 
at a radial distance of 150 Rj. The results are shown in Fig. 9. One must be 

careful not to construe this figure as being what the spacecraft actually ob- 

served. Such should be the case only for arcs of slight curvature, i.e., 
short time duration, otherwise Jovian rotation and lo revolution effects need 
be taken into account. In Figure 9 only the model arcs of small curvature can 

be directly compared to the data as explained above. This includes the lo 

I 

longitudes from approximately 100" to 230" when the spacecraft is in the re- 
spective system III longitude range from approximately 60" to 300®. At both 
high and low lo longitudes the modelled arc signatures change. For some arcs 
the time duration becomes greater than thirty minutes. For others the curva- 
ture becomes quite high. These morphological changes may be due to azimuthal 
assymmetries in the magnetic field 0-4 model. The low curvature arcs occur 
when the angle, 4>, at Jupiter between the spacecraft and lo differs only 
slightly from 90® {by say no more than 20"). 

We find that our model of low curvature arcs with » 90° agrees 
reasonably wall with the data set— the V2 examples cited in Figure 8 for days 
79-182 and 79-197 are typical cases. Low curvature arcs as observed by 
Voyager 1 on day 79-63 from about 950 to 1030 are shown in Figure 2, At 1005 
in this Figure the spacecraft longitude is given as 152" and the lo longitude 
at the same time is 238". Referring to Figure 9 we find good agreement 
between the observed arcs and the model-produced low-curvature arcs at the 
same approximate spacecraft longitude and lo longitude. 
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was assumed that the spacecraft was at 150 Rj and at a latitude of 3.2**. 

Again, each point on each model curve represents the spacecraft longitude at 
which the chosen frequency intersects the approximate spacecraft latitude at 
the time of observation. Hence each point was determined by intersecting the 
approximate spacecraft latitude with the model-generated emission cone surface 
at a radial distance of 150 Rj, !lhe results are shown in Fig. 9. One must be 

careful not to construe this figure as being what the spacecraft actually ob- 

served, Such should be the case only for arcs of slight curvature, i.e., 
short time duration, otherwise Jovian rotation and lo revolution effects need 
be taken into account. In Figure 9 only the model arcs of small curvature can 

be directly compared to the data as explained above. This includes the lo 

longitudes from approximately 100” to 230” when the spacecraft is in the re- 
spective system III longitude range from approximately 60® to 300®. At both 
high and low lo longitudes the modelled arc signatures change. For some arcs 
the time duration becomes greater than thirty minutes. For others the curva- 
ture becomes quite high. These morphological changes may be due to azimuthal 
assymmetries in the magnetic field 0-4 model. The low curvature arcs occur 
when the angle, (}), at Jupiter between the spacecraft and lo differs only 
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slightly from 90° (by say no more than 20®). 

We find that our model of low curvature arcs with ()j « 90® agrees 
reasonably well with the data set — tlie V2 examples cited in Figure 8 for days 
79-182 and 79-197 are typical cases. Low curvature arcs as observed by 
Voyager 1 on day 79-63 from about 950 to 1030 are shown in Figure 2. At 1005 
in this Figure the spacecraft longitude is given as 152® and the lo longitude 
at the same time is 238® . Referring to Figure 9 we find good agreement 
between the observed arcs and the model-produced low-curvature arcs at the 
same approximate spacecraft longitude and lo longitude. 
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OF POOR QUALFaY 


Our emission mechanism assumes that DAM is the result of Doppler-shifted 
gyro-emission from accelerated, field-aligned electrons. As noted by Omidi 
and Gurnett (1982) the condition that must obtain at relativistic gyro- 
emission resonance (condition for maximum wave growth rate) is 


P 



0 ) 0 ) 

n cos (i) n cos (b 



( 1 ) 


where "n" is the phase index of refraction and the subscript "g" refers to the 
gyrofrequency. Solving (1) for dcjj/df yields 

d4. -0iq(2 - 

df 4Tm p sin (|) f^ 


This expression suggests that the observed peak in 4i(f) must be explained geo- 
metrically since no physical values of the variable parameters will allow 
d(l*/df > 0. 

For both cases in Figure 10 the lo flux tube intersects Jupiter at a lon- 
gitude where little azimuthal assymmetry in the 0-4 B-field model is present 
(particularly so for the VI example at lo III longitude «180® [during arc 
emission the lo longitude changes by a few degrees] ), The results seem to be 
best explained by noting that since the source positions for the highest fre- 
quencies along the lo flux tube are closer to the planet by about 1 Rj, in- 
creased refraction of the higher frequencies refracts the emission cones so 
that the spacecraft observes the higher frequencies "inside" the emission 
cones of the lower frequencies in a three dimensional ray path. 
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Further analysis of Equation 1 has allowed us to suggest the most 
probable value of the Doppler shifting present for the higher-curvature arcs. 
Goldreich and Lynden-Bell give the ratio f/fg for the proposed Doppler-shifted 
gyro-emission as 


- . 1 . 1 ^ 


(3) 


For a reasonable range of values of p, Equation (1) allowed us to set limits 
on the range of cjj. Table I contains calculated values of the maximum for a 
given p. 


TABLE I 


p 

f/fg 

^’max 

0.05 

1 .0025 

85.7° 

0.1 

1 .01 

81 .4° 

0.15 

1 .023 

76.95 

0.2 

1 .04 

73.2° 

0.25 

1 .0625 

71 .9° 

0.30 

1 .09 

65.6° 


Table I indicates the relativistic beaming effect and strong dependence 
of 4i on the ratio, f/fg/ allowed by the resonance condition, equation (1). 

Recall that the low curvature arcs in most cases can be modelled ade- 
quately for (|> w 90° wconstant. The higher-curvature arcs, on the other hand, 
require smaller wave normal angles by about 15° -20° for typical examples. 
Hashimoto and Goldstein (1982) have also found from their ray-tracing results 
that as f/fg increases, smaller initial values of 4> (w70°) gave better fits to 
the data. This suggests that the higher-curvature arcs are produced by a more 
energetic electron field-aligned beam. Our results indicate, for the cases 


15 


studied in detai.l, that 4^n,ax ~ indicating that « .16 or f/fg 1.02. 

The low curvature arcs require little or no beaming to explain their flat 

(low-curvature) signatures. While we can say little about the low curvature 

arc electron energies, the higher-curvature arc electron energies can be found 

using p = .16 to be about 6.5 keV. 

As in the case of the low curvature arcs, some of the higher-curvature 
arcs, cannot be explained by the model. Occasionally in the data set a series 
of arcs of higher curvature occur when the angle between the spacecraft and 
lo, <J) 180“. For thei^e cases (one good example is VI day 79-34 at 12 - 16 

hrs. [Pig. 12] ) it is not possible to ray trace with the present 0-4 mag- 
netic field model, because an angle of (j> » 180“ would require the ray path to 
pass through Jupiter in order to be observed by the spacecraft. One possible 
way out of this dilemma is the introduction of a new B-field model such as the 
Smith et. al. (1975) D-4 model which might allow such ray paths. This is the 
case because the D-4 model has the geographic pole of Jupiter completely out- 
side the lo-auroral oval, so that field lines with a foot on one hemisphere 
(for instance the north) might pass over the pole and above the opposite hemi- 
sphere. This would allow a source point located on one side of Jupiter at the 
foot of an lo flux tube to propagate to the opposite side of Jupiter from lo 
( (j) « 180“). 
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SUMMARY I 

We have presented Jovian DAM ray-tracing results using the 0-4 magnetic 
field model and a Jovian empirical-fit plasma model that includes the lo 
torus. The assumed emission mechanism requires that: (1) the radiation is 

R-X mode; (2) the source region is at the foot of an lo flux tube and the 
source mechanism is Doppler-shifted gyro-emission, and (3) the emission cone 
is hollow. The ray-tracing code is three-dimensional and based on the Stix 
cold-plasma formulation of the index of refraction and Haselgrove's set of 
first order differential equations. 

We have shovm that DAM ray paths refracted by the ambient plasma and 
local B-field can produce arc-like signatures as observed from the VI and V2 
spacecraft explained by intersecting the emission cones with the spacecraft 
trajectory. Because of the dependence of f^-x plasma and B-field parame- 
ters, a dispersion of frequencies occurs which duplicates the observed signa- 
tures in the data, in many, but not all, cases studied. 

Our model arcs explain the vertex-early and vertex-late arcs of low cur- 
vature as due to very small Doppler-shifted gyro-emission emitted initially at 
constant wave normal angle, (|j » 90®, at all frequencies. There is good quali- 
tative agreement for most of the low curvature arcs examined when the space- 
craft is relatively far from Jupiter to minimize spacecraft motion effects. 
Some data does not match the model; however, as in cases where the model sug- 
gests no arcs should exist because the spacecraf t-Io angle, <|), is much larger 
than 90°. For most examples of low-curvature arcs (cf. nested arcs of large 
frequency extent) the model successfully predicts the general shape and occur- 
rence of the data examined. Decametric arcs of higher curvature seem to rep- 
resent a different class of event. 
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In the data aet one frequently observes either isolated or nested groups 
of arcs that have a larger curvature that cannot be explained if one assumes 
» constant m 90®. We have found one possible mechanism to adequately model 
this class of signature is to allow the wave normal to vary with frequency. 
This result yields a 4>(f) relationship that seems to be typical of a large 
number of higher-curvature arcs. The tj>(f) curve that typifies this class of 
arc is peaked around 10 MHz and falls off on either side with almost constant 
but different slopes the positive slope being somewhat larger than the nega- 
tive slope. A semi -quantitative analysis of the (|i{f) relationship based on 
the resonance condition (1), allows us to use the observed maximum value of c|) 
for the modeled higher-curvature arcs, tj^max “76®, to suggest p « 0.16. This 
value, based on the Doppler-shift expression (3) implies the maximum energy of 
higher-curvature arc-producing electrons is »6.5 keV, a reasonable value com- 
pared to suggested potential differences of several hundreds of kV across lo, 
a possible source of electron acceleration along the field lines (cf, Shawhan, 


1976 ). 
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CONCLUSIONS I 

We conclude from our results that many qualitative features of the DAM 
arcs are explained by assuming a gyro-emission mechanism along an lo flux 
tube. Resultant emissions are refracted by the ambient plasma and B-field to 
produce the arc-like signatures observed by Voyagers I and II. Both the shape 
and location of the data arcs examined are predicted. Our results imply that 
higher-curvature arcs may be due to more energetic electrons whose ultimate 
gyro-emission is more relativistically beamed to smaller wave normal angles. 
There are still unexplained observations. Our model has not yet explained 
some low-curvature emissions that occur for 4> significantly different from 90" 
and some higher-curvature emissions for which ^ » 180"(Io independent arcs). 

It may be necessary to incorporate a different B-field model such as the Smith 
et. al, D-4 model to explain the latter result, or invoke an electron source 
other than lo. The Warwick concept of B-field anomalies and/ or the Boischot 
model of arcs of the same frequency and source position producing an interfer- 
ence pattern cannot be ruled out for some arc events . We hope to resolve the 
discrepancies by future improvements of the model. 
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WHISTLER MODE RAY TRACING 

The purpose of this study is to determine the maximum extent in latitude 
and in longitude of the Jovian high latitude lightning source region. Prom 
this, one can obtain a lower limit to the Jovian surface lightning density: 


lightning density 


no. of strokes 
area 


The ray tracing code is the same as described in the previous section — a 
3-dimensional code based on the Haselgrove (1955) ray tracing equations and 
using the cold plasma index of refraction given by Stix (1962). The magnetic 
field model is that of Acuna and Ness (1976) with coefficients set to repre- 
sent a dipole with a magnetic moment of 4.225 gauss Rj^. The model of the 
magnetospheric plasma density is again the Sentman-Goertz plasma model 
( Sentman-Goertz , 1977). The lo plasma torus model was obtained by a spline 
fit interpolation of the density contours published by Wairwick et al. (1979). 
This model is symmetric about the magnetic equator with a maximum density at 
the center of the torus of about 2500 cm"3 (Menietti fi Gurnett, 1980). The 
total plasma density model used in the ray-tracing computer code is shown in 
Figure 3. The Jovian ionosphere used in this model is an empirical fit to the 
Voyager I ingress data published by Gehrels (1979) (Fig. 13). The maximum 
ionospheric density occurs at an altitude of approximately 2000 km. Prom the 
ray tracing code with the empirical ionospheric model included, the index of 
refraction at h « 2000 km was found to be 

"max 6.08 . 

Using Snell's law let us consider an atmospheric lightning burst refracting 


into the ionosphere and obtain limiting conditions 
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Consider two points A and A' in the Jovian atmosphere which represent the 
maximum latitudinal extent of a lightning burst (Fig. 14). The distance 
between points A and A' represents the maximum extent of the lightning source 
region. For the limiting case consider a ray from point A to the south and 
tangent to the bottom of the ionosphere. Assuming the atmospheric index of 
refraction to be 1 r we find, from Snell’s law, 

sin A f* 1 

where nj[^ = the ionospheric index of refraction 

A = angle between the wave normal vector and the normal to the surface 
at some point of the ray path (see Fig. 14). 

For the limiting case, nj[ = %ax 

^ ~ t 

so 

sin Ac = 1/6.08 
A = 9.47“ . 

So rays from the atmosphere reaching the bottom of the ionosphere (C and C' ) 
with A > Ac = 9.47" will not reach the magnetosphere, but will be internally 
reflected. In a grossly exaggerated fashion Fig. 14 shows rays "1" and ”2” 
from points A and A' respectively which both trace to point D in the lo torus. 
By ray tracing, we have found the points A and A' which will allow the ray 
patlis of rays 1 and 2 to intersect at D in the Jo torus. 

In this study rays at a constant frequency of 3 kHz were launched from 
high Jovian latitude ( «60“ ) . For the latitudinal study the rays were launched 
in a meridian plane with the angle A (from Figure 14) fixed at 9.47“. It was 
assumed the satellite was located in the lo torus on the magnetic equator at r 
= 6.0Rj. By trial and error it was found that rays launched from latitude = 
58.6", A = 9.47“, cross the magnetic equator at r = 6.02 Rj (corresponding to 
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ray 2 in Figure 14). Rays launched from latitude - 62.6“, A ■ 9.47“, cross 
the magnetic equator at r « 5.99 Rj (corresponding to ray 1 in Figure 14). 

The horizontal extent of points A and B is 4.0“ and corresponds to a surface 
distance of about 4978 km. This is a large distance but represents the maxi- 
mum extent detectable since we have taken the limiting case. 

The longitudinal extent of the lightning source was investigated in a 
similar fashion. The ray source point was moved in longitude from a given 
reference lightning source longitude. Rays were launched at A 9.47“ while 
varying angle 4», ^e angle between the wave normal vectox and the magnetic 
field vector and varying p, the angle the wave normal makes relative to the 

meridian plane. By varying both the ray launch longitude and angle p while 
keeping A ^ 9.47“ ray paths were found that crossed the magnetic equator at 
r w 6.01 Rj. For larger excursions in longitude the ray path will not inter- 
sect the satellite position within the torus. So it is concluded that the 
maximum extent of the lightning source region in the longitudinal sense at 
this latitude is ~1.9“, corresponding to a horizontal distance of 2360 km. At 
latitude = 62.7“ (the top of the latitudinal extent), longitude = ± .85“ and p 
= 113“, rays cross the magnetic equator at r = 5.999 Rj. For larger excur- 
sions in longitude at this latitude the ray path will not intersect the satel- 
lite position within the torus. For this latitude the maximum longitudinal 
extent of the lightning source region is 1.7“ corresponding to a horizontal 
distance of 1991 km. To reiterate, these are the maximum extents in longitude 
because we have considered the limiting case. If one crudely assumes a rect- 
angle with maximum latitudinal extent of 4980 km and maximum longitudinal ex- 
tent of 2400 km, the maximum area of detectable lightning source on the Jovian 
surface at high latitude is 1.2 x 10^ km2. Hence a lower limit to the light- 
ning source density can be found if one knows the number of lightning strokes 
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occurring within this area, A summary of the latitude and longitude excursion 
computer runs is given in Table II, 


—A,., 


Source 

A 

TABLE II 

Latitude Excursion Runs 

Values Magnetic Equator Crossing 

Latitude Longitude Radius Longitude 

26.4® 

VD 

0 

0 

-9.466® 

o 

. 

CO 

in 

0® 

6.01 5Rj 

0® 

5.06® 

90® 

+9.466® 

62.6® 

0® 

5.987RJ 

0® 




Longitudinal Excursion 

Runs 




Source 

Values 


Magnetic Equator Crossing 



P 

A 

Latitude 

Longitude 

Radius 

Longitude 

26.4® 

94® 

-9.41® 

58.6® 

0.95® 

6.009RJ 

0® 

5.0® 

113® 

9.86® 

62.7® 

0.89® 

5.9985RJ 

0® 


In order to evaluate the wave energy dispersion, four rays were launched, 
each with a constant wave normal angle, ()) (either (j> * 1®, (|< » 5®, or (}> « io®). 
Each ray was launched at a different angle p ( p = 0® , 90®, 180®, 170®), The 
area intercepted by this envelope of rays at the top of the ionosphere and at 
the equator was measured. The rays' source points were the center of the 
latitudinal and longitudinal lightning source region determined above. The 
index of refraction and area at both points (a and b in Fig. 15) as well as 
the ratio of these two areas, A2/A1 are listed in Table III, 
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jL 

Index of Refraction 

Area (Hm^) 

A 2 A 1 

A 


1.76 

Al «• 6,7 X 10^ 




3,9 X 103 

B 

10 

47 < n < 53 

A2 * 2.6 X 10^ 


A 

5® 

1.76 

Al » 1.7 X 1 q5 

3.8 X 103 

B 

5® 

32 < n < 82 

A2 “ 6,4 X 10® 


A 

10® 

1.77 

Al - 7.3 X 10® 

3.8 X 10® 

B 

10® 

51 < n < 96 

A2 “ 2.8 X 10® 





SU14M)VRY II 



The 

whistler mode ray-tracing studies have demonstrated that wj assuming 

a source 

point at high Jovian latitude and ray-tracing 

through the 

magnetosphere to the lo torus, an 

estimate of the lightning source extent on 

the 

Jovian surface can be made. This estimate is much 

smaller than that made 


by Lepping et, al (1981) who have performed the same study by starting at the 
lo torus and ray-tracing back toward the planeto This latter study achieved a 
surface area at least three times higher than ours. The discrepancy,, we feel, 
is due to the high degree of wave refraction near the Jovian surface and 
consequent Landau damping effects which were not considered by Lepping et al. 
These damping effects ,nre not as important for rays launched from the surface 
due to the higher phase velocity and smaller wave normal angle of the waves at 
the same point near the surface on the outbound trace compared to the inbound 
trace. 

Our results can now be utilized to estimate the lightning source rate at 
Jupiter and, from the energy dispersion calculated, determine the approximate 
Jovian radio power output due to lightning bursts on the surface. 
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NEW 10 TORUS MODEL 

The lo toruB contours presented by Warwick et. al. (1981) were a result 
of a first loo); at the Voyager 1 data. To draw tJie contours the assumption of 
torus symmetry with respect to the Jovigraphic equator was made. As a firtit 
approximation this is not a bad assumption; however, Hill et. al. (1974) have 
shown that the centrifugal equator is a much better assumed central axis of 
the lo torus. More recently, and post-Warwic); lo torus model, Bagenal and 
Sullivan (1980) have presented a much more sophisticated loo); at the lo torus 
including the addition of ion contributions and both hot and cold species. 

The results of their wor); are shown in Figure 16 which depicts the lo torus 
electron density contours* This figure, when compared to Figure 1, immedi- 
ately shows the effect of the hot electron species which are depicted as lobes 
or "horns" at approximately 5 Rj. Since the Bagenal-Sullivan model is a more 
accurate description of the torus it has been introduced into our new ray- 
tracing code. 

To reproduce the results of Bagenal and Sullivan linear interpolation of 
table values are made numerically in two dimensions: along the radial direc- 

tion and along the magnetic field line. The table values present electron 
density along the magnetic field line as a function of L-shell. A grid of 
points is represented between radial distances of 5 and 8,5 Rj and up to 1,5 
Rj above the centrifugal equator at 0,1 Rj intervals. The points in the grid 
have been put on the meridional plane where the magnetic and rotational equa- 
tors are aligned (km « 291®) so that the plasma is symmetrically distributed 
about these equators. Azimuthal symmetry is also assumed. 

An Euler transformation allows us to transform from Jovigraphic coordi- 
nates to Jovian magnetic coordintes which are tilted with respect to the 
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latter by the complement of the angle, a, shovm in Figure 17. It is a three- 
dimensional transformation, so that any longitude can be chosen. The angle y 
in Figure 17 can be shown for a dipole field to be 

Y « tan-"* {2/3 tan a/(1 + Jl + 8/9 tan^ a)^) (1) 

The arc length along the magnetic field line can be shown to be 

® “ ■'*1 TT^ li * ji * ■ i’-" It* 

where u *» cos0. 

The numerical number density can then be found by evaluating 5n/6r and 
5n/d0, where the derivatives can be evaluated directly from table values after 
first noting that 

5n/S0 = (5n/5s) (ds/50) 
and Ss/50 is evaluated from (2) above. 

SUMMARY III 

The Bagenal-Sullivan torus has been incorporated into the new ray-tracing 
code, and will provide a much better basis for the ray-tracing of low 


frequency radiation in the lo torus. 
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Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 
Figure 6 


Figure 7 


Figure 8 


Figure 9 


Figure 1 




FIGURE CAPTIONS 

• A large isolated arc double -valued in frequency is seen beginning 
at about 1155 UT at both high and low frequencies with a vertex at 
about 1235 UT. 

. Frequency versus time spectrogram for VI -79-63 depicting a series 
of low curvature nested arcs especially clear between about 950 UT 
and 1030 UT. 

. Model Jovian Magnetosphere. Contours of Plasma frequency are from 
Sentman and Goertz (1978), B-field from Acuna, and Ness (1976), lo 
torus from Warwick, et al (1979). 

■ Calculated Radiation Sheets for Two Northern Hemisphere Sources: 

2 MHz and 20 MHz, and for Two Southern Hemisphere Sources: 2 MHz 

and 10 MHz. The radial shell at six Jupiter radii (6Rj) is shown 
in mercator projection. The source points are designated by "X" 
and all rays are launched at 90“ with respect to the B-field. The 
lo torus is shown in longitudinal cross section, shaded. 

, Calculated Emission Cones for the Same Source Points as in Pig. 4., 
but now Intersecting the lOORj Radial Shell. 

. Calculated Emission Cones for Wave Normal Angles (WNA's) of 60“ , 

80“ , 100®, and 120“. The mercator projection is taken at lOORj, 
top panel: 2MHz, bottom panel: 25 MHz. 

I Calculated Emission Cones Emanating from Sources Located at 
Different Longitudes. All rays have been launched at 90“ , wave 
frequency 25 MHz, mercator projection at lOORj. 

, Frequency versus system III '65 Longitude. The points designated by 
squares are model points. The numbers above the model curves 
depict the approximate lo flux tube longitude chosen for the 
source points. The event on the left is V2-79-197; the right 
event is V2-79-182. 

, Frequency versus system III '65 Longitude showing the model arcs for 
each 10“ interval of longitude. The lo longitude is labelled on 
each arc. The wave normal angle at the source point for all these 
arcs was 90“ . 

I. High Curvature Arcs, Double-Valued in Frequency, are shown here for 
Vl-79-71 (left) and V2-79-188. These arcs are fairly typical of 
high curvature arcs and were chosen because of their isolated 
nature making them easier to measure. 
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Figure 


Figure 


Figure 


Figure 


Figure 


Figure 


Figure 


11. Wave Normal TVngle, versus frequency found empirically by varying 
<|) until the mo'^el arcs matched the observed arcs for VI — 79-72 
(top) and V2-79-188 (bottom). 

12. Frequency versus time spectrogram for VI -79-34. Several arcs occur 
between 1347 UT and 1527 UT while the angle, (Ji, between the space- 
craft and lo is approximately 180“. 

13. Altitude versus Electron Number Density for the Jovian Ionosphere 
as Determined during Voyager 1 Ingress and Published by 
Gehrels (1979). 

14. A Schematic Diagram of Whistler Ray Paths "1" and "2" from Points 
A and A' Respectively which Both Trace to Point D in the lo torus. 
The points A and A' represent the maximum latitudinal extent of a 
lightning burst on the surface of Jupiter. "N" is the normal to 
the surface at A and A' . C and C’ designate the bottom of the 
ionosphere; A is the angle between the wave normal vector and the 
normal to the surface at some point of the ray path. 

15. A Schematic Diagram Depicting the Ray "Cone" Produced by Four 
Whistler Waves Stimulated by Lightning in the Center of the Jovian 
Lightning Source Region (as defined in the text) . The dispersion 
of this ray cone is crudely found by comparing the area A, at 
point "a", the top of the ionosphere, to the area A 2 at point "b", 
the magnetic equator. 

16. Density Contours of the lo Torus as Published by Bagenal and 
Sullivan (1981). The vertical axis is a measure of the distance 
above the centrifugal equator (as defined in the text) . The 
dashed lines depict the inbound and outbound trajectories of VI . 

17. A Schematic Diagram showing the Jovicentric Angular Velocity, Si, 
tilted with respect to the magnetic moment vector. 111, by the 
compliment of the angle a. The angle y is between the centrifugal 
equator and the magnetic equator. The arc length, s, is along the 
(dotted) dipole magnetic field line between some reference point at 
angle 6 and the magnetic equator. The actual field lines (solid) 
are shown distorted by the centrifugal plasma pressure. 


30 


ACKNOWLEDGEMENTS 

We would like to thank the GSFC radio astronomy group for supplying us 
Voyager I and II PRA data. We would also like to thank Jerome De La Cruz 
his assistance in programming. This work was supported by NASA contract 
NAS7-100 through Cal Tech/JPL Contract 956026. 


with 

for 


3 


ORIGINAL PAGE f9 
OF POOR quality 



VWfiGER 1 79/ ,63 MflR 4 6: 0: 0 TO 13: 0: 0 

n?COlJENCY TItt SPeCTROGRftM POWER FUJX FILTD?- 0.0U 


C »Ci». ' I 1 

Cf QUALITY 









1 ^ Iff 

5HK 


'▼ ... *ir. 


• a 

0 UM m ^ 

»•■ ■ ^b"fSr- 

'•■V * .w . A ’ P^.C J?.a \<W - M ^ 

• iZ ,£l 

. ■ -*s W •- 

>tX)T*7<L 


S»8S 

- !2Ki8?g 


ir-” 4.1.-“ rf's 


?■ ^^2^ 
5- P§i§ 


in in ^ 

^ (\J N 


Q VP 9 h- 
- ^ w ^ in 

N ^ 

in (i) «-* «H 

M fo o) ^ 

” KD *-« 


» ,« F-> I 1.3 

»- CL 














LATITUDE LATITUDE 




f=2MHz 




FREQUENCY= 26MHz SHELL = 100R, V=90 





SYSTEM III (1965) LONGITUDE 


o 



SYSTEM III LONGITUDE 



sc R=150 RJ LAT=3.200 FG/F RATIO= 1.005 PSI ==90 





OF POOR QUAUTY 


O 

T" 

CO 



JL 

CJ 


CM O 00 CO ^ CM O 

CM CM r- -I- T- t- r- 


(ZHI/'J) AON3n03Ud 


00 


o 

CM 


(O 




CM 


O 

CO 

CM 


SYSTEM III LONGITUDE 



■ Wi^ER 1 75^ 34 FEB 3 11:30: 0 TO 17: 0: 0 

^TCiXCNCY TI(C SPECTROGRAM POWER FLUX FILTER- 0.01X 


ORIGIN/Vl page i3 
OF POOR QUALITY 





cm 





ORIGiNAL 
OF POOR QUAU 



t*. 


magnetic equator 


HEIGHT (RJ). 





RADIAL DISTANCE (RJ) 




ORIC-' 

OF pc 







